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ABSTRACT 

Using the unique X-ray imaging capabilities of the Chandra observatory, a 
2006 observation of Cygnus X-3 has provided insight into a singular feature asso- 
ciated with this well-known microquasar. This extended emission, located ~ 16" 
from Cygnus X-3, varies in flux and orbital phase (shifted by 0.56 in phase) with 
Cygnus X-3, acting like a celestial X-ray "mirror". The feature's spectrum, flux 
and time variations allow us to determine the location, size, density, and mass of 
the scatterer. We find that the scatterer is a Bok globule located along our line 
of sight, and discuss its relationship to Cygnus X-3. This is the first time such a 
feature has been identified with the Chandra X-ray Observatory. 

Subject headings: X-rays: binaries — X-rays: individual (Cygnus X-3) — X-rays: 
ISM 



Introduction 



At a distance of 9 kpc (IPredehl et al. 1120001). Cygnus X-3 is an unusual microquasar in 
which a Wolf-Rayet companion (Ivan Kerkwijk et al. Ill992l ) orbits a compact object with an 
orbita l period of 4.8 hours. It is a strong radio source routinely producing radio flares of 1 to ~ 
20 Jy (IWaltman et al. II1995I ). It has been shown to produce radi o jets and the radio emission 



correlates with both the soft X-ray and hard X-ray emissions flMioduszewski et al. 112001 
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Miller- Jones et al. 112004 ; ISzostek et al. \\200A McCollough et al. Ill999f ). In 2000 Chandra 



obse rvations found exten ded X-ray emission that is believed to be associated with Cygnus 
X-3 (jHeindl et al. II2003I ). In this paper we analyze a 2006 Chandra o bservations of Cygnu s 



X-3 and reexamine the previous Chandra observations discussed in iHeindl et al. I (120031 ). 



In particular we take a careful look at the timing and spectral properties of the extended 
feature and how they related to Cygnus X-3. 



2. Observations and Analysis 



Between 1999-2006 Cygnus X-3 was observed six times by Chandra using the ACIS- 
S/HETG with a Timed Event (TE) mode and covered a range of source activity (see Table 
1). These observations provided grating spectra of Cygnus X-3 and on-axis zero order images 
which allow for spatial and spectral analysis of Cygnus X-3 and its surroundings. 

The primary observation used for this analysis was a 50 ksec quenched state observation 
(OBSID: 6601 hereafter referred to as QS) obtained during a period of high X-ray activity. At 
the time of the observation the RXTE/ASM (2-12 keV) count rates were ~ 25 — 30 cts s -1 , 
Swift/BAT (15-50 keV) had an average count rate of ~ 0.0 cts s _1 , and the Ryle radio 
telescope (15 G Hz) showed radio fluxes of ~ 1 mJy, all va l ues typical of a Cygn us X-3 



quenched state (ISzostek et al. M2008I : iMcCollough et al. 111999 : IWaltman et al. 111996 



J. This 



observation has the longest duration of any of the observations and the feature has the 
greatest number of photons of any of the observations. 

Since this analysis involves a region near a bright point source we also need to determine 
the impacts of pileup on the data. This can be determ ined by looking at the average number 
of counts per detection island (9 pixels region: see the I Chandra POG I (120111 )) per observing 
frame. We do this by taking a series of annuli centered on Cygnus X-3 for the highest count 
rate observation (QS). Each annulus has a radial thickness of 2 ACIS pixels (~ 1") with 
the readout streak regions excluded. We then sum the counts in each annulus and divide 
by the number of observation framed and the of number detection islands (the area of the 
annulus in pixels divided by 9 pixels). In Figure 1 is a plot of counts/frame as a function 
of radial distance from Cygnus X-3, with the solid line represents the entire observation 



1 The number of observing frames can be determined either by extracting the keyword TIMEDEL from the 
event file header and dividing the exposure time by its value or using the acisf06601_000N001_statl.fits file 
found in the Chandra secondary data products which contains a record of all of the frames in the observation. 
This file must also have applied all the time corrections and filtering that were applied to the event file and 
correct CCD selection needs to be made. 
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and the dotted line the times of peak count rate. From iDavis I (120071 ) we have taken the 
counts/frame values for which one would expected pileup of 1% , 5%, and 10% and plotted 
them in Figure 1. We see that at radial distance of greater than ~ 8" pileup should not be 
an issue and will not impact our analysis. 

In our analysis of these observations, we used version 4.3 of the CIAO tools. The 
Chandra data retrieved from the archive were processed with ASCDS version 7.7.6 or higher. 
All event files were filter by their good times intervals and barycenter corrected using the 
CIAO tool axbary. 



2.1. The Feature and its Characteristics 



For each observation, we have extracted a zeroth-order image between 1 — 8 keV binned 
at the nominal ACIS-S resolution (~ 0.49"). Figure 2 and 5 shows images for the QS ob- 
servation. In each individual image, there is a bright, unresolve d core with strong r adial 
point-spread function (PSF) wings and a strong scat t ering halo ( IPredehl et al. 1 120001 ). In 
each observation the feature reported by lHeindl et al. I ( 120031 ) is present at the same location. 
An analysis of the QS observation shows the feature (RA : 20 h 32 m 27.1 s , DEC : +40°57'33.8") 
lies at an angular distance of 15.6" from Cygnus X-3 at an angle of 68.5° from the orien- 
tation of the observed radio jets of Cygnus X-3. The feature is extended and was fit using 
CIAO/SHERPA with a 2D-Gaussian with axes of 3.6" and 5.5" and a position angle of 78.7°, 
measured counterclockwise with the top of the image being 0°. If the feature and Cygnus 
X-3 are at the same distance, then their separation is 2.2 Dg light-years (D 9 : distance in 
units of 9 kpc). This feature is present in all observations of Cygnus X-3 done throughout 
the Chandra mission (1999 to present). 



2.2. Temporal Behavior 



The feature is located amidst a high background region. This high background is due 
to the telescope PSF and a dust scattering halo. The majority of the photons in the back- 
ground have energies above 2 keV and as such arise from scatte ring off micro-roughness in 
the telescope optics; see the discussion in IPredehl et al. I ( 120001 ) . In order to examine the 



temporal variability of the feature we need to subtract the background flux from region of 
feature. To do this we used segments of annuli centered on Cygnus X-3 for the feature and 
background regions to extract the light curve for the feature and the background (see Figure 
2). The annuli were chosen to maximize the number of counts for the feature as well as 
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have the feature and background regions sample the same radial region of the PSF/dust 
halo while avoiding the readout streaks. All data extraction was done with the CIAO tool 
dmextract. 



2.2.1. Phase Folded Analysis 

The phase-folded (1-8 keV) light curves of Cygnus X-3, the background region, and 
background subtracted feature region are shown in Figure 3. As might be expected the 
background which is mainly due to the scattered PSF emission demonstrates the same 4.79 
hr orbital variation with the slow rise and rapid observed drop in Cygnus X-3. This is 
also reinforced in Figure 4 (top panel) which show no observed lag in the cross-correlation 
between Cygnus X-3 and the background. However, the feature surprisingly shows the same 
orbital modulation but with a phase shift of ~ 0. 6. Phase-selected images reveal how this 
feature varies relative to Cygnus X-3 (see Figure 5), which can be seen more dramatically in 
the movie which accompanies this paper (see online version). This is also confirmed in the 
cross-correlation between Cygnus X-3 and the feature seen in Figure 4 (lower panel) which 
show a clear lag that peaks around 9460 seconds which corresponds to a phase lag of ~ 0.55. 

Fitting Phase Folded Light Curves: To address questions of a possible issues with back- 
ground subtraction (over-subtraction) we took the light curve from feature annulus and fit 
it using the following model: 



Cif[i] = C bkgd [i] + a * C bkgd [shift (i,j)]. (1) 

Where Cif[i] is the count rate for the feature's region in phase bin i, Cbkgd [i] is the count 
rate from background region (scaled to the feature's region size) in phase bin i, and the last 
term a* Cbkgd [shift (i, j)] is the count rate of the background region scaled by a a nd shifted by 



j in p hase. The fit parameters are a and j. We used the IDL routine MPFIT (lMarkwardt 
20081 ) to fit the light curve. We did multiple fits of the light curve using a number of different 
phase binnings and found good fits for all binnings with consistent fit values (see Table 2). 
Our best fit a gave a = 0.29 ± 0.01 and phase shift(j) = 0.56 ± 0.02 (see Figure 6). 

As a further test we replaced the second term in Equation 1 with a constant term which 
was used as a fit parameter. We found no acceptable fits for any of the binnings (see Table 
2). This result is consistent with the feature varying with the same period as Cygnus X-3 
but shifted in phase. 



Phase Image: Finally, as another way to check to see if background subtraction could 
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be an issue for each Chandra detected event, a "Cygnus X-3" phase value was determined 
from its arrival time. Photons falling into certain phase ranges were broken in separate 
"phase" images. These images were assigned a certain color and combined to form a color 
coded phase image. The bands were: (red) 0.3-0.63, (green) 0.63-0.96, and (blue) 
. 96-0 . 3. The image is shown in Figure 7, note the blue color of the feature. This indicates 
that bulk of the photons are arriving in the 0.96-0.3 phase range. If the feature was constant 
we would expect it to be white (equal amounts of each color). It is also important to note 
that no background subtraction was done in this approach and hence there is no issue with 
the background subtraction creating a false time/phase variation of the feature. 



2.2.2. Longer timescale Variability 

To examine the temporal behavior of the feature on longer timescales we determined 
the fiujfl of Cygnus X-3 and the feature for each of the Chandra/ ACIS observations (see 
Table 1). When Cygnus X-3's flux is plotted versus the flux from the feature there appears 
to be a correlation (see Figure 8) where the dashed line is a linear fit to the data. A 
Pearson's correlation test of the data yields a correlation coefficient of 0.98 indicating a 
linear correlation between the feature and Cygnus X-3. 



2.3. Spectrum 



The annuli used to extract the light curves (see Figure 2) were also used to extract 
spectra for the feature and background region. The final background-subtracted spectra 
of the feature contained from 78 to 3216 counts in the 1-8 keV band. Table 3 shows fits 
to the spectrum using simple absorbed power law and blackbody models all of which yield 
acceptable fits to the spectrum. 

All of the Chandra/ ACIS observations were HETG grating observations. The ± first or- 
der HEG spectra were combined and fit for each observation (see notes i n Table 1). The con- 
tinuum for Cygnus X-3 is complex and dependent on the state of activity (jHjarmarsdotter et al. 
20081 . l2009l ; iKolionen et al. Iboioh . In the X-ray (0.5-10 keV) the continuum can be mod eled 
by partial covered disk blackbody during flaring/quenched states (IKoljonen et al. 1120101 ) and 
during the quiescence/transition states we found that continuum was best approximated by 
an absorbed power law. In all cases a large number of spectral features were added (see 



2 The spectral fits and corresponding fluxes estimates are given and discussed in section 2.3 and 3.1. 
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Table 4) to improve the spectral fit. 

We note that all of the spectra of the feature are more absorbed and very steep/soft, rel- 
ative to the corresponding Cygnus X-3 spectra (see Table 3 & 4). As Cygnus X-3 transitions 
from a quiescent (hard state) to flaring/quenched state (soft state) its spectrum becomes 
softer. The spectrum of the feature is shown to follow suit and becomes steeper/softer as 
Cygnus X-3's spectrum does. 



3. Nature of the Feature 

To understand the feature's nature the following must be explained: (a) it is clearly 
extended; (b) its flux varies in phase with Cygnus X-3 with a phase shift of 0.56; (c) the flux 
from the feature shows a correlation with the flux from Cygnus X-3; (d) the phase-averaged 
flux of the feature is ~ 10~ 3 of Cygnus X-3's flux; (e) the time variation is 4.8 hrs but the 
separation between the feature and Cygnus X-3 is at least 2.2 D 9 light-years; and (f ) its 
spectrum is heavily absorbed and lacks hard X-ray flux relative to Cygnus X-3. 

Jet Emission: It is natural to try to associate this feature with the Cygnus X-3 jet 
emission. However, Cygnus X-3 was in a quenched state throughout the Chandra observa- 
tion and for several days before and after, during which jet activity is strongly suppressed. 
Additionally, over a three ye ar period prior to this observations Ryle/AMI-LA observed no 



radio flare exceeding 0.5 Jy (jPooley Il201ll ). Furthermore, in earlier Chandra observations 
where there is jet activity, the feature is fainter. If the X-ray emission of the feature was 
due to the jet one would expect this to be synchrotron emission and one would not generally 
expect the spectrum to be so steep and heavily absorbed. Also the misalignment of the 
feature relative to the jets observed in the radio complicates this picture. 

Jet Impact Are a: Problems with this being a jet impact area have been noted by 



Heindl et al. I (J2003J) (location relative to the radio jets, jet precession, and jet collimation). 
In addition the observed flux correlation between the feature and Cygnus X-3 provides prob- 
lems given a likely separation of 2.2 Dg light years. The strong phase modulation of the 
feature (varying the same as Cygnus X-3 by a factor of two), combined with a periodicity 
exactly matching that of Cygnus X-3, would be difficult to understand. One would expect 
the continuing impact of the jet would give rise to a brighter constant flux from the feature 
and drastically reduce the modulation which is observed. 



Wind Interaction: The strongest arguments which were made by iHeindl et al. I (120031 ) 
for the feature's nature is that it is due to a wind/ISM interaction. The feature's distance, 
flux correlation and phase modulation with Cygnus X-3 make such a model difficult to 
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reconcile with the observations. In this interpretation the feature is created over timescales 
that are long (~ 2000 years) compared to the observed orbital modulation. Also the direct 
change in the flux of the feature with Cygnus X-3 flux is difficult to reconcile with such a 
model. 



3.1. Scattering Solution 



A natural explanation for the feature's time variable behavior is that it is a result 
of scattering from a cloud (which acts as a kind of interstellar X-ray "mirror") between 
Cygnus X-3 and the observer. This explanation would naturally lead to the observed phase 
difference between light curves as a difference in the path length of the scattered photons 
( iTrumpler fc Schonf elder 1 1 1 9 73h . This could also explain the flux correlation (see Figure 
8) and overall flux difference (~ 10~ 3 ) since only a small fraction of the total flux will 
be scattered toward the observer. In figure 9 a diagram of the expected geometry for the 
scattering from a cloud is shown. 

The spectrum of the feature can also be modeled as a scattered version of Cygnus X-3's 
spect rum. At these energies scattering will modify the spectrum by A * E~ 2 (ISmith et al. 
20021 ) . due to the reduced scattering efficiency at higher energies (lOverbeck Ill965l ). There 
will also be an additional reduction at low energies caused by absorption in the cloud (N c i) 
and multiple scattering. We can then model the feature's spectrum as due to scattering as 



Slf = e - CT ( E > NH W*A*E- Q *S cont . (2) 

Where is the scattering model for the feature, e _(r ^ NH ^ is the additional absorp- 
tion due to the cloud (modeled using phabs in XSPEC), A * E~ a represents the high energy 
attenuation due to scattering (modeled using plabs from XSPEC), and S CO nt is Cygnus X-3 
continuum model determined from the grating data (see Table 4). The resulting fit parame- 
ters for the observations are shown in Table 5 and the fit to the QS observation's spectrum 
is shown in Figure 10. The resultant fluxes (1 — 8 keV) for the feature and Cygnus X-3 are 
given in Tables 4 & 5. The flux ratio of the feature to Cygnus X-3 are all in the range of 
(4 - 6) x 10" 4 . 

This gives rise to a natural explanation of the observed spectral differences. The loss of 
the flux below 2 keV is the result of additional absorption in the cloud and the reduced flux 
at higher energy simply reflects the drop in scattering efficiency. 

Distance to and Size of the Feature: Assuming that this feature is due to an individual 
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cloud the time delay is determined by the distan ce to the cloud (jTrumpler fc Schonfelder 



19731 ; iPredehl fc Klose 111996c iPredehl et al. Il2000l ). The resulting time delay can be written 



as: 

At = = !. 15 e^. (3) 

2c 1 - x 1 - x w 

Where At is the time delay (in seconds), D is the distance (in kpc) to the source, is observed 
angular distance (in arcsec) from the source, x is the fractional distance of the scatter to the 
observer (see Figure 9), and c is the speed of light. From the observed phase observed phase 
offset, we know that time delay is given by At = (0.56 + n)t cx3 where t cx3 = 17.25 ksec is 
the observed orbital period of Cygnus X-3, but we have an ambiguity in the total number 
orbital period offsets (n). Using At = 0.56 *t cx3 = 9.66 ksec , D = 9 kpc, and 6 = 15.6", 
the resulting fractional distance is x = 0.79. This means that the cloud is close to Cygnus 
X-3 (within 1.9 kpc) and if n > 0, this distance could be less. 

This degeneracy (n) can be removed using the fact that the feature is extended. From 
equation (3) we would expect that the delay scattered photons experience increases as a 
function of angle from Cygnus X-3. If we take the inner and outer edges of the feature to 
be 13.8" and 17.4" respectively then for x = 0.79 we would get a time delay of ~ 1.2 hrs 
across the feature. For locations closer to Cygnus X-3 we would expect the delay across the 
feature to increase. To test this we have taken the extraction and background annuli (see 
Figure 2) and divided into an inner and outer set of annuli (11.2" to 15.6" and 15.6" to 20.0" 
respectively). In Figure 11 we show a cross correlation of the inner and outer light curves 
(for the QS observation with 5 minute time resolution) which shows a significant lag (at the 
99 % confidence level) at 0.9 hrs. This would correspond to the n = case. In Figure 12 
are the phase folded light curves of the inner and outer region in which one can see that the 
light curve for the outer region is lagging the inner by ~ 0.2 in phase. This corresponds to 
a lag of ~ 1 hr. This puts the feature at a distance of 1.9 kpc from Cygnus X-3 making the 
observed dimensions of the feature 0.12 by 0.19 parsecs. 

Density of the Cloud: From the spectral fit to the feature for the QS observations 
we have an absorbing column density of the feature of 5.0^' 7 x 10 22 cm -2 . This value is 
consistent with column density determined from spectral fits to the other observations (see 
Table 5). The column density of the feature can also be estimated from the flux ratio of the 
feature to Cygnus X-3 using the following relationship (see Appendix A for derivation) 



^ = N H (lf) 

cx3 



yrtano!! tana 2 cos 2 ( 
fl-x) 2 
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x 



K r S^e-^Hwf™ a" 3 ' 5 ( d <^,9 a) j ^ 



.i=g,si 



(4) 



The flux ratio is equal to the product of three quantities. The first (Nn(lf)) is the 
column density of the feature, the second is a solid angle term which relates to what fraction 
of Cygnus X-3's flux is intercepted by the feature, and final term is a scattering term which 
is a measure of the flux scattered by the dust in the cloud (this is solved for by numerical 
integration over the energy range and grain distribution). The last two terms are depend on 
x, the fractional distance between the observer and Cygnus X-3. Figure 13 shows a plot of 
the last two terms and their product as function of x. The solid line is for scattering due to 
silicates and the dotted line for scattering due to graphite. Equation 4 can be solved for the 
column density of the feature as a function of x as is shown in Figure 14. The parameters used 
to create these curves are given in Table 6. Also included are lines representing the column 
density from the spectral fit of QS with its uncertainties and a vertical line representing its 
location from the observed time delay. What we find is good agreement with the column 
density determined from the spectral fits and the column density (Ny = 3.6 x 10 22 cm -2 ) 
necessary to produce the observed flux ratio. 

If we assume that path-length along the line-of-sight through the cloud is similar to the 
observed dimensions of the feature [(3.7 — 5.9) x 10 17 cm] we can make an estimate of the 
density of the feature. Taking the column density to be 5.0l 2 'y x 10 22 cm -2 we arrive at a 
density range of (0.6 — 1.9) x 10 5 cm -3 making the feature a dense molecular cloud. 

Mass of the Cloud: From the estimate of the density of the cloud and the size determined 
from the X-ray measurements we can make an estimate of the mass of the cloud. Using the 
simplifying assumption of a spherical cloud with a diameter of between 0.12 — 0.19 pc with 
a density of 10 5 cm -3 we arrive, from X-ray observations alone, at an estimate of the mass 
of the cloud to be 2 - 24 M . 

From th e cloud's size, density, and mass the featu re has all of the characteristics of a 



Bok Globule (IBok fc Reilly 1119471 ; IClemens et al. Ill99ll ). but instead of seeing this as a dark 



obscuring feature in the optical we see it shining in scattered X-rays. 



4. Relationship to Cygnus X-3 

What is the relationship of this feature to Cygnus X-3? Three possibilities present 
themselves: 

(1) Random Alignment: Cygnus X-3 and the feature both lie in the Galactic plane 
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(lii = 7 9.845°, bjj = 0.700°) . T he Cygnus X region that hosts Cygnus X-3 is rich in molecular 



clouds ( ISchneider et al. 1120061 ). So this may be just a chance alignment. If so this gives us 
insight into the nature and structure of molecular clouds in the ISM. In this case we would 
be looking across three Galactic arms (with Cygnus X near the Local Spur, the feature in 
the Perseus arm at ~ 5 kpc, and Cygnus X-3 in the Outer Arm at 7 — 9 kpc). Bringing 
Cygnus X-3 closer, to 7 kpc, does not greatly change the distance estimate to the feature or 
the need for three star forming regions along the line of sight. 

(2) Supergia nt Bubble Shell: These structures have been observed in other galaxies 



( IKim et al. Ill999l ). They have typical radii of 0.5-1.0 kpc and are driven by the radiation 



and outflows from OB associations, supernovae and their remnants, and XR Bs. Molecular 
clouds have also been found to exist in these shells (jYamaguchi et al. Il200ll ). Cygnus X-3 
is a high-mass X-ray binary and likely still resides in such an OB association. This gives a 
natural explanation of the feature's location along our line-of-sight. Given the distance of 
the feature from Cygnus X-3 this w ould be a bubble comparable to the HI shell found in 
NGC 6822 fide Blok &: Walter I boooh . 



(3) Microquasar Jet-Inflat ed Bubble: Within NGC 7793 a powerful microquasar is 
driving a 300 pc jet-inflated bubble ( IPakull et al. Il2010|). Cygnus X-3 is a microquasar whose 



jets a ppear to be aimed along our line of sight f lMioduszewski et al. 11200 It iMiller- Jones et al. 
20041 ) . It is possible that rapid cooling near the working surface of the jet, in the shell of 
the cocoon, may allow a dense molecular cloud to form. This would explain the nature 
of the feature as well as its alignment with Cygnus X-3. Although it should be noted 
that there is research which sugge st that the jet may not be close to the line of sight 
( jMartf. Paredes fc Peracaula 1120011 ). which would make this a less likely option. 



Finally it should be noted that a combination of both (2) and (3) may be possible. 
The radiation and outflows from the OB association may create a large low density cavity 
in which the microquasar jet can more easily propagate. Evidence for large-scale cavities 
surrounding other microquasars has been noted (lHao fc Zhang 1 l2009h . This could explain the 
large distance of the feature from Cygnus X-3 (1.9 kpc) and reduce the energetics necessary 
to produce it. The feature may be located at the place where jet interacts which the wall of 
the cavity. 



5. Conclusion 



This feature and its temporal relationship to Cygnus X-3 have unveiled the unique 
interaction between a microquasar and its environment. It has given us a tool to probe the 
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nature and structure of molecular clouds, providing information on their size and shape, 
possibly due to the microquasar interaction or the presence of ordered magnetic fields in the 
ISM. It has also given us our first X-ray view of a Bok Globule. 

To date this is the first such feature found with Chandra. If the feature is indeed due to 
a microquasar interacting with its environment then we would expect there to be very few. 
This would be due to the limits of small angle scattering in the X-ray and the need for the 
lobes and associated molecular clouds to be aligned close to our line of sight. Depending on 
the nature of these sources the best candidates would be high mass X-ray binaries (because 
of their young age and hence likely relationship with an OB association and star forming 
regions) with a relatively short orbital periods. 
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A. Derivation of Scattering Relations for Cygnus X-3's and the Feature 

It is possible to derive some of the scatter's properties by comparison of the fluxes of 
Cygnus X-3 and the feature . This derivation is similar to that done for the scattering halo 



intensity of ISmith fc Dwek I ((1998). The geometry being used can be found in Figure 9. If 
we take the unabsorbed luminosity of Cygnus X-3 as a function of energy to be L cx3 (E) then 
the X-ray luminosity at the feature is given by 

ME) = . (A1) 

Where cr(E) is the X-ray absorption cross section, Nn(r s ) is column density along the 
path (r s ) between Cygnus X-3 and the cloud. 

The photons scattered into a solid angle dQ by a single scatter for a source of luminosity 
L S (E) is given by 
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P(E)=L s (E)^dfi. 



(A2) 



Where da K /dfl is t he di fferential scattering cross section (see iMathis fc Lee I ( I199ll ); 



Mauche fc Gorenstein I (|l986h ). 



For a telescope with collecting area A' we can chose df2 such that A' = r Q 2 dfi, where r Q 
is the distance from the scatter to the observer. 

The photon count rate that the observer will detect from scattering from a single dust 
particle as 



C s (E)^L s (E)e-^M(^)^ 



L cx3 (E)e-( E )[ N H(-)+ N H(r )] 



47rr 2 V dfi / r n 2 



Where Nn(r ) is the column density between the scatter and the observer. 

The feature has been fit with an elliptical Gaussian with semi-major and semi-minor 
axes of ri and r 2 respectively. For a distance of xD one can use the angular measurements 
of the semi-major and semi-minor axes, a\ and 02 respectively, to give the surface area of 
the feature as Ay = 7rrir 2 — 7rx 2 D 2 tan«i tan«2- For a scattering region thickness of l s the 
scattering volume of the feature is given by = Ajfl s . For a dust grain number density of 
n g we can determine the total count rate and flux for the feature as 

v ( v\ C ^(E) n g (a)L cx3 (E)e-( E )[ N H( r s) + N H(ro )] / d(7 x A lf l s 

In general n g will be a function of dust radius a. For our case we will consider n g to be 
uniform spatially at the location of the feature. Integrating over the dust size distribution 
gives the flux for the feature as 

Fi,(E)= X w UJ^ da - (A5) 

We can simplify the above equation by noting that the observed angle 9 for the feature 
is very small and the overall angular dimensions ( ot\ and 0:2) of the feature are small. In 
this case the path traveled by the scattered photon will be very close to the path traveled by 
the unscattered photon. Because of this the total hydrogen column density for both paths 
should be the same except for the additional column density of NH(lf) along the scattered 
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path due to the feature. If we take Nh(cx3) to be the column density between the observer 
and Cygnus X-3 then the column density along the path of the scattered photon can be 
written as 



N H (r s ) + N H (r G ) = N H (cx3) + N H (lf) . (A6) 
F cx3 (E) the observed absorbed flux from Cygnus X-3 can be written as 

F«s(E) = L -^^ • (A7) 

Finally from the scattering geometry (see Figure 9) we have 

1 cos( 



D( 1 — x) 



(Ah 



Where x is the projected distance along path between the observer and Cygnus X-3. 
The measure angle of the feature (9) and the scattering angle (9 S ) are related to the projected 
distance x by 9 = (1 — x)9 s . Using the above relations we arrive a the following substitution^] 

L cx3 (E)e-( E >[ N H(r S )+NHM] F cx3 (E)e-( E ) N HWcos 2 (g s -g) 

4vrr s 2 ~ (1-x) 2 ' 1 j 

If we integrate over energy bandpass and replace F cx3 (E) by F cx3 S(E) where F cx3 repre- 
sents the measured flux from Cygnus X-3 and S(E) is its spectral form (normalized to one) 
we can write the flux relationship of the feature to Cygnus X-3 as 



Fif Ti tan«! tan a 2 l s cos 2 

F^ = (1-x) 2 

S(E)e" 



E2 /"a m a> 

<r(K)N H (If) 



dcr s (E, a, 9 S ) 
i).,(rt) I y — I da<lL . 



(AlO) 



Assuming a MRN grain size distribution ( iMathis. Rumpl. fc Nordsieck l ll977t IWeingartner fc Draine 



20011 ) then we have 



3 We note that this equation is similar to Equation A5 in lSmith fc Dwek~ ( 1998h . However, in that paper, 
one factor of cos(6> s — 9) was omitted in error, although this makes no difference to either their or our final 
results. 
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n„ a 



J2 n j 



-3.5 



(All) 



Where nh is the hydrogen number density of the cloud, a is the radius of the grain and 
NJj is are the normalization in (grains/H atom)///m for graphite (g) and silicates (si). If we 
also note that n h l s is simple the column density of the cloud N H (lf). Substituting n g (a) in 
A10 gives us 



cx3 



N H (lf) 



7r tan ai tan a 2 cos 



x 



5>; 



.l=g,Sl 



E 2 

d , S(E)e 

Ei 



(1-x)* 

-a(E)N H (lf) 



-3.5 



dcr s (E, a, 9 S 



dadE 



• (A12) 
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Fig. 1. — Plotted are the counts/frame for the QS observation as a function of radial distance 
from Cygnus X-3. The solid line is for the entire observation and the dotted line is for the 
times of the highest count rates (Cyg nus X-3 phas es 0.6-0.7). The counts/frame for various 
levels of expected pileup are given (IDavis 1 120071 ). The long dashed vertical lines is the 
location of the feature and the shaded region give the feature's radial extent. Pileup should 
not be an issue for the observations used in this analysis. 



-18 - 




Fig. 2. — Segments of annuli centered on Cygnus X-3 were used for the feature and back- 
ground regions to extract the light curve and spectrum for the feature (blue) and the back- 
ground (green). The inner edge of the annuli is 11.2" from Cygnus X-3 and the outer edge 
is 20.0". 
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Fig. 3. — Phase folded 1-8 keV light curves of Cygnus X-3, the background, and the feature 
from the QS observation. Top: Phased fold light curve of Cygnus X-3. The data were taken 
from the ± first order HEG and MEG spectra. Middle: Phase folded light curve of the 
background region. Bottom: Phase folded light curve of the feature from QS (background 
subtracted). It can clearly be seen that the feature exhibits the same slow rise and rapid 
drop that one sees in Cygnus X-3 but with a phase lag of ~ 0.6. 
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Time Log (sec) 

Fig. 4. — Cross correlation plots: The dashed vertical line corresponds to a lag of half of 

Cygnus X-3's period and the dot-dashed corresponds to a lag of Cygnus X-3's period. Top: 
The cross correlation between Cygnus X-3 ( grating data) and the background region (1-8 
keV) using 600 second time samples. Note there is no lag between them. Bottom: The cross 
correlation between Cygnus X-3 ( grating data) and the feature (background subtracted) 
using 172 second time samples (~ 0.01 of Cygnus X-3's period). Note the lag of 9460 
seconds which corresponds to a phase shift of ~ 0.55. 



Fig. 5. — These phase selected images were created from the zero-order image Chandra 
grating observation (QS) using photons with energies between 1-8 kev. The count range and 
color scales are identical for both images. The bright readout streak caused by the ACIS 
CCD readout and the "cratering" of the central source due to pileup are both instrumental 
effects, left: Phase range 0.96-0.3 image. The feature is prominent, right: Phase range 
0.5-0.8 image. The feature is very weak. Online electronic edition: To better visualize 
the phase relationship between the feature and Cygnus X-3 a movie was created. For all 
of the events detected in QS a Cygnus X-3 orbital phase based on their arrival time was 
determined. The data was filtered to only include 1-8 keV energy photons. For the image 
a 256 by 256 pixels region centered on Cygnus X-3 was used. A set of images were created 
based on phase intervals, starting at 0.0 phase and with a duration of 0.20 of Cygnus X-3's 
phase. Each image was successively shifted by a 0.01 of Cygnus X-3's phase until images 
for the full Cygnus X-3 period were created and compiled to form the movie. The duration 
of the observation (50 ksec) corresponds to 2.96 orbits of Cygnus X-3, so phase coverage is 
relatively uniform. The duration of each image (0.2 phase) corresponds to an integration 
time of ~ 10 4 seconds. 
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Fig. 6. — Fits to the phase folded 1-8 keV light curves of the feature with no background 
subtraction. Top: Phase folded light curve of Cygnus X-3 (50 phase bins) with a fit to the 
light curve (dotted line). Middle: Data minus model fit divided by the data error bars. With 
± 1 a error bars. Bottom: The model of the phase folded light curve of the feature from the 
fit. 
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Fig. 7. — A color coded phase image of Cygnus X-3 region. Note the blue color of the 
feature. This indicates that bulk of the photons are arriving in the 0.96-0.3 phase range. It 
is also important to note that no background subtraction was done and hence there is no 
issue with the background subtraction creating a false time/phase variation of the feature. 
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Fig. 8. — Plotted is the 1 — 8 keV flux of Cygnus X-3 (determined from the grating data) 
versus the 1 — 8 keV flux of the feature for each of the ACIS observations using the scattering 
model. The dashed line is a linear fit to the data. 
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Fig. 9.— Scattering Diagram Definitions: Above is shown a diagram of how the 
scattering is taking place with various parameters labeled. X : Location of the X-ray source 
(Cygnus X-3). S : Location of the feature. O : Location of the observer. D : Distance from 
the observer to Cygnus X-3. x : Distance to the scatter along the path to Cygnus X-3. 9 : 
Observed angle of the feature from Cygnus X-3. 9 S : The scattering angle of the X-rays from 
Cygnus X-3 by the feature. r s : Distance from Cygnus X-3 to the feature. r Q : Distance 
from the feature to the observer, h : Distance from the the feature to Cygnus X-3's line of 
sight to the observer. 
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Fig. 10. — The extracted spectrum of the feature, taken from QS, is shown above. Also 
given is the fit (in red) and residuals of a scattering model. A good agreement is found 
between the data and model. 
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Fig. 11. — A plot of the cross correlation between the inner (11.2" to 15.6") and outer (15.6" 
to 20.0") regions of the feature (both background subtracted). The time resolution was 5 
minutes for the light curves. The dashed horizontal lines corresponds to a 99% confidence 
level. The dash, dot-dash, dots-dash, and large dash vertical lines correspond to lags one 
would observe across the feature for n = 0, 1, 2, and 3 respectively. Note the prominent 
peak at ~ 1 hour corresponding to n = lag. 
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Fig. 12. — A plot of the phase folded light curves of the inner (solid: blue online) and 
outer (dotted: red online) regions of the feature (both background subtracted). There is a 
noticeable lag in outer relative to the inner of ~ 0.2 in phase which corresponds to ~ 1 hr 
in time. 
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Fig. 13. — Plots of last two terms and their product of the flux ratio (see equation 4) as a 
function of x. Top: Plot of the solid angle term which is a measure of the flux the feature 
intercepts. Middle: Plot is of the scattering terms which takes into account what fraction 
of the X-ray flux is scattered to the observer. The solid lines is for silicates scatters and the 
dotted line is for graphite scatters. Bottom: This final plot is the product of the solid angle 
term with the sum of the two scattering terms. 
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Fig. 14. — A plot of the scattering column density necessary to produce the observed fea- 
ture/Cygnus X-3 flux ratio (see equation 4) vs x. The vertical line (large dashes) is location 
for the feature determined from the time delay. The horizontal dashed lines represent the 
scattering column density determined from spectral fit and the shaded region between the 
dashed-dotted lines represent it uncertainties. We find good agreement between values found 
for the flux ratio and the spectra fit. 



Table 1. Cygnus X-3 HETG Chandra Observations 



ObsID 


Date (MJD) 


Instrument 


Data Mode 


Exp Mode 


Exposure (ksec) 


Count Rate a (cts/sec) 


State 


101 b 


51471 


ACIS-S 


FAINT 


TE 


1.95 


35.8 


t/mf 


1456 (obi 0) b 


51471 


ACIS-S 


FAINT 


TE 


12.12 


27.7 


t/mf 


1456 (obi 2) 


51531 


ACIS-S 


FAINT 


TE 


8.42 


18.0 


q/qi 


425 c 


51638 


ACIS-S 


GRADED 


TE 


18.54 


88.2 


fs / mrf 


426 c 


51640 


ACIS-S 


GRADED 


TE 


15.68 


71.0 


fi/mrf 


6601 


53761 


ACIS-S 


FAINT 


TE 


49.56 


106.6 


h/qu 



Note. — The states are given as k/s where: fcare t hose of Kolionen et al. I ( 2010l) [q: Quiescent, t: Transition, fh: FHXR, fi: 
FIM, fs: FSXR, and h: Hypersoft] and s are those of Szostek et al. ( 2008h [qi: Quiescent, mf: Minor flaring, su: Suppressed, 
qu: quenched, mrf: Major flaring, and pf: Post flaring]. 

a This is calculated by taking the total number of 1-8 keV events in the observation and dividing it by the exposure. 

b The observation end time of obsid 101 is identical to the observation start time of 1456 (obi 0). For the purpose of the 
analysis in this paper these two obsids will taken to be a single observation. 

c Obsids 425 and 426 were done in alternating exposure mode. The shorter frame time observation (0.3 sec) resulted in an 
additional 0.4 ksec (425) and 0.3 ksec (426) exposure. For this analysis only the observations with the longer frame time (1.8 
sec) were used. 
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Table 2. Phase Folded Light Curve Fitting Parameters 



Model 


Binning 


Amplitude 


Phase Shift 


X 2 /dof 




10 


0.29 ±0.01 


0.60 ±0.10 


(9.24/8) 




20 


0.29 ±0.01 


0.55 ±0.05 


(23.74/18) 


Phase Shift 


40 


0.29 ±0.01 


0.55 ±0.03 


(50.82/38) 




50 


0.29 ±0.01 


0.56 ±0.02 


(51.17/48) 




100 


0.29 ±0.01 


0.56 ±0.02 


(104.60/98) 




10 


0.069 ± 0.002 




(66.34/9) 




20 


0.068 ±0.002 




(80.81/19) 


Constant 


40 


0.068 ±0.002 




(108.07/39) 




50 


0.068 ±0.002 




(110.38/49) 




100 


0.067 ±0.002 




(166.39/99) 



Note. — The Phase shifted model is the one give in Equation 1. 
The Constant model is the same as the Phase shifted model except 
that the shift term has been replaced with a constant which is used 
as a fit parameter. 
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Table 3. Spectral Fit Parameters for the Feature 



Obsid: 


101+1456(0) 


1456(2) 


425 


426 


6601 


NetCounts a : 


163 


78 


750 


697 


3216 


Absorbed Power Law: 


N h (10 22 cm- 2 ) 

r 

Flux b 


7 q+16.7 
'•°-7.3 

1 7+ 4 - 2 

L - 1 -2.9 

6.6 x 10- 13 


6.0l 36 5 6 

9 S+ 4 - 5 

3.6 x 10- 13 


12.5^; 2 7 
4.8l 2 ;° 
2.6 x 10~ 12 


13.3+^ 
2.6 x 10- 12 


10.6l 2 - 8 

4 n+0.6 
^• u -0.6 

3.7 x 10- 12 


X 2 /dof 


5.8/8 


7.7/13 


14.5/13 


17.2/13 


15.2/17 


Absorbed Blackbody: 


N h (10 22 cm' 2 ) 
T (keV) 
Flux b 


1 o+13.3 

o 4 +o.9 
7.1 x 10- 13 


9 0+24.4 
z -<J-2.2 

1 1 +L6 
1 - 1 -0.6 

3.8 x 10- 13 


7 q+4.5 
' -°-2.7 

7 +0 - 2 

u - ' -0.2 

2.5 x 10- 12 


i-o±8:i 

2.6 x 10- 12 


5.611; 3 

0.918:1 

3.5 x 10- 12 


X 2 /dof 


5.9/9 


7.3/13 


13.6/13 


16.4/13 


15.5/17 



a Net number of counts in 1-8 keV part of the spectrum. 
b Measured flux in ergs sec -1 cm -2 in the 1-8 keV band. 
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Table 4. Continuum Spectral Fit Parameters for Cygnus X-3 



Obsid: 



101+1456(0) 1456(2) 



425 



426 



6601 



N 



10 22 cm" 2 ) 

r 

T (keV) a 
pc (10 
f, 



:!2 cm 2 ) 



pc 

Flux b 

X 2 /dof 



4 12 +0 - 17 



0.72 



-0.05 
-0.05 



1.69 x 10~ 9 
4151.55/8660 



6.18^ 

o.63i8:« 



1.11 x 10" 9 
1668.70/4319 



3.41 



-0.25 
-0.26 



2.25 



-0.04 
-0.04 



o qi+0.28 
°- y± -0.24 
-0.06 
-0.06 



0.74^ 
5.86 x 10~ y 
3269/4317 



4 34 +0 - 47 



9 99+0.05 
z - zz -0.04 

4-58±gJg 

0-851^ 
4.95 x 10~ 9 
2502.15/4317 



9 7 r+0.11 
z - ' d -0.10 



2 00 +0 02 

z - uu -0.02 
o 99+0.22 
°- zz -0.19 

62 +0 - 04 

u - uz -0.04 

6.13 x 10~ 9 
6016.95/4317 



Note. - The Chandra grating spectral are rich in spectral features (jPaerels et al. Il2000l ). 
To achieve acceptable fits we included a large number of spectral features (emission lines, 
absorption lines, edges, and radiative recombination continuum). 

a Temperature of a disk blackbody. 

b Measured flux in ergs sec -1 cm -2 in the 1-8 keV band. 



Table 5. Spectral Fit Parameters for Scattering Model for the Feature 



Obsid: 


101+1456(0) 


1456(2) 


425 


426 


6601 


N d (10 22 cm" 2 ) 


r 4+ 15 - 5 

°-^-4.19 


o o+ 7 - 8 




4.01ft 


5.01 2 -? 


A 


1 4+0-5 

-L-^— 0.5 


1 5+°- 7 


1 l +0 - 2 
i - ± -0.2 


i-i±8:l 




a 


2.0* 


2.0* 


2.0* 


1.811J 


i q+o.6 

0.6 


Flux a 


6.0 x 10~ 13 


3.9 x 10~ 13 


2.7 x 10~ 12 


2.6 x 10~ 12 


3.6 x 10~ 12 


X 2 /dof 


6.4/8 


7.8/13 


15.8/13 


16.7/12 


14.7/16 



* Fixed. 

a Measured flux in ergs sec -1 cm~ 2 in the 1-8 keV band. 



Table 6. Scattering Flux Calculation Parameters 



Parameter Value 



F cx3 (1-8 keV) 9.12 x 1CT 1 photons chi'V 1 (6.13 x 1(T 9 ergs chi^V 1 ) 

(1-8 keV) 5.6 x 10~ 4 photons cm _2 s _1 (3.6 x 10~ 12 ergs cm" 2 s _1 ) 

e 15.6" 

a x 2.77" (5 Sigma Full Width: 5.54") 

a 2 1.82" (5 Sigma Full Width: 3.65") 

D 9 kpc 

p {graphite) 01 2.24 g cm -3 

p (silicate) 0, 3.5 g cm -3 

N d (graphite)* 7.41 x 1(T 16 

N d (silicate) a 7.76 x lO" 16 H g a r t a Q i °% 

a min a 0.005 fi 

0.25 p 

Ex 1.0 keV 

E 2 8.0 keV 



The dust grain parameters were taken from 



Weingartner &: Draine 



(12001( 1 . 



